Electronic states around a vortex core in high-T c superconductors 

based on the t-J model 
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Electronic states around vortex cores in high-T c superconductors are studied using the two- 
dimensional t-J model in order to treat the d-wave superconductivity with short coherence length 
and the antiferromagnetic (AF) instability within the same framework. We focus on the disappear- 
ance of the large zero-energy peak in the local density of states observed in high-T c superconductors. 
When the system is near the optimum doping, we find that the local AF correlation develops inside 
the vortex cores. However, the detailed doping dependence calculations confirm that the experimen- 
tally observed reduction of the zero-energy peak is more reasonably attributed to the smallness of 
the core size rather than to the AF correlation developed inside the core. The correlation between 
the spatial dependence of the core states and the core radius is discussed. 



There are several experimental results in high-T c su- 
perconductors which have not been explained in the 
conventional BCS eL„2_ a 2-wave superconductivity. The 
electronic structure inside the vortex core is one of 
the most interesting issues among them. The conven- 
tional theory for d-wave vortices based on Bogoliubov- 
de Gennes (BdG) mean-field theory predicts a large and 
broad peak at the Fermi energy in the local density 
of states (LDOSL so-called zero-energy peak (ZEP), at 
the vortex core|l|. However, scanning tunneling spec- 
troscopy (STS) spectrum in one of the high-T c materi- 
als, BSCCO, giving directly the LDOS around the vortex 
core, shows only a small-double peak structure at ener- 
gies ±7meV0. A similar situation was also observed in 
YBCO compounds^. 

To resolve this discrepancy, there have been proposed 
several theoretical attempts, that is, a d x 2_ y 2 +s stateQ, 
a d x 2_ y 2+idx^\ state^H, an antiferromagnetic (AF) vor- 
tex core [3, la la], a staggered flux state|lij, and a vortex 
core with small /cf£o EH E3 where kp is the Fermi wave 
number and £o is the coherence length. These theories 
are, in greater or lesser degree, based on strongly corre- 
lation effects. Among them, we consider the AF vortex 
core and the effects of small k-p£o (the small vortex core) 
because they have some experimental grounds mentioned 
below. 

The AF vortex core was first predicted using a phe- 
nomenological Ginzburg-Landau theoryQ assuming the 
SO(5) symmetry 1 13|. Later this possibility was con- 
firmed in a microscopic calculations li| In the t- 
J model, it was shown that the AF core is stabi- 
lized near the doping rate d ~ 0.1, where the sys- 
tem is close to the AF instability^). Simultaneously 
it was observed that the electron density approaches 
half- filling inside the AF core@. After these calcula- 
tions, there are several experimental studies suggesting 
the existence of the enhancement of AF correlations or 



AF vortex core, such as neutron scattering experiments 
in LSCO compounds near optimal doping [lit Hfl IT^ . 
and NMRdllU 

experiments in underdoped or 
near-optimally doped YBCO. Therefore it seems natural 
that the AF vortex core is also realized in the BSCCO 
sample 2] which is optimum doped or slightly overdoped, 
and that the absence of the large ZEP is due to the mean- 
field AF gap (or Mott gap) developed inside the core. 
Actually some theoretical studies for the AF vortex core 
have predicted the LDOS without the ZEPH@, or at 
least, large-double peak structure [ill |22| . 

The double peak structure is also found in the small 
vortex cores (small /c_f£o) by u sing the BdG equations 
for continuum d-wave model|lll|l2j. In fact, the STM/S 
data has revealed that a radius for the vortex core in 
BSCCO is on the order of 1(M 2] , which is consistent with 
the short coherence length. Therefore, the small vortex 
core is also a possible explanation for the absence of the 
ZEP, but the lattice effects, which play an important role 
when kp^o is small, need to be clarified. 

Motivated by these backgrounds, in this paper we per- 
form a detailed study on the electronic states around vor- 
tex cores in high-T c superconductors, taking account of 
the discrete lattice effects, smallness of the vortex core, 
and the AF instability. To investigate their effects within 
the same framework, we use the t-J model which has a 
phase diagram consistent with experiments^.^ and en- 
ables us to study the reasonable doping dependence. We 
find that the reduction of ZEP is more reasonably at- 
tributed to the smallness of the core size rather than to 
the AF correlation developed inside the core. It is re- 
markable that a simple-minded picture of AF core does 
not play important roles in the suppression of LDOS. 

The Hamiltonian of the t-J model is written as 
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(1) sites including two vortices. In this case, by the appro- 



(ij) 



in the standard notation where means the sum- 

mation over nearest-neighbor pairs. The Gutzwiller's 
projection operator Pq is defined as Pg = IIi(l — 
ni\riii). The uniform magnetic field is introduced in 
terms of Peierls phase of the hopping term as ty = 

t exp [f- c J?' A • dr) 0, 0- The BdG equation based on 
the extended Gutzwiller approximation is 
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where i + t represents the nearest neighbor sites of the 
site i, and a = ±1. The renormalized parameters tff , 
J<f and /i?f are determined from cluster calculations [24l 
which reproduce the variational Monte Carlo results 
They have the following forms: 
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where g t (i, j) , g x s y '{%, j) and g z s {i,j) depend on the lo- 
cal expectation values Ay = ^ ( c lt c i4, ~~ c ll c Jt^' ^ = 
{c\ a Cj a ), and 771$ = jj(nj| — n^) 25]. For example, 
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where £y (my) is an average of hole density (AF mo- 
ment) for the corresponding bond and ay is a factor close 
to one which depends on Ay , Xij > $ij an d % [24| • These 
somewhat complicated renormalized parameters are nec- 
essary for treating the AF and d-wave superconductivity 
in the same framework. If one uses simple slave-boson 
mean- field theories |26j|. the AF order parameter is over- 
estimated and thus one cannot discuss the AF vortex 
core because the AF order exists even in the bulk. In the 
present method, AF order does not exist away from the 
core for S > 0.10, which is consistent with the numerical 
results 0. 

To make the Peierls phase compatible with vortex lat- 
tice symmetry, we need a magnetic unit cell with 2N x N 



priate choice of gaugeyj,^], the order parameter Ay has 
a translational symmetry with respect to the magnetic 
unit cell. Here we take N =22-26. We solve numeri- 
cally the BdG equation and carry out an iteration un- 
til Ay , Ay , Xij i Xij > n icr are determined self-consistcntly. 
We take J/t = 0.3 and examine the various values of the 
doping rates S = 0.10-0.15, that is, from near optimum 
doped to overdoped region. 




FIG. 1: (a) Spatial dependence of the amplitude of d-wave 
(upper) and s-wave (lower) components of the superconduct- 
ing order parameter around the vortex core. The doping rate 
is S = 0.15. (b) The hole density n%. (c) The LDOS obtained 
at the vortex core center (the thick line) and obtained without 
the magnetic field (the thin line). 

First, let us look at the results for the overdoped re- 
gion. Figure ^ shows the superconducting order param- 
eters, the hole density, and the LDOS obtained around 
the vortex core for 6 = 0.15. We can see the common 
features found in several theoretical studies within the 
weak-coupling d-wave superconductivity 0, that is, fairly 
large vortex cores, slightly induced s-wave component|^|, 
the increase of the hole density inside the core, and the 
ZEP in the LDOS. We note that the local AF order inside 
the core does not exist in this case. 

The situation drastically changes when the system ap- 
proaches the optimum doping. Figure 2 shows the results 
for 5 = 0.11 in the same manner as in Fig. ^ Similar re- 
sults for S = 0.10 have been shown before 8], but here 
much larger system-size is achieved. From Fig. [21 we can 
see that the size of the vortex core is fairly small, that is, 
the core having a radius of about 3 lattice spacings is re- 
alized. The magnitude of the induced s-wave component 
is also small as in the case of S = 0.15. Figure|2tb) shows 
the spatial dependences of the local AF moment \rrii\ and 
of the hole density n\. It is found that AF correlation 
develops and the hole density decreases inside the vortex 
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FIG. 2: (a) Spatial dependence of the amplitude of d-wave 
(upper) and s-wave (lower) components of the superconduct- 
ing order parameter around the vortex core. The doping rate 
is 5 — 0.11. (b) The hole density rv£ (upper) and local an- 
tifcrromagnetic moment (lower), (c) The LDOS obtained at 
the vortex core center (the thick line) and obtained without 
the magnetic field (the thin line). 



core, in sharp contrast with the S — 0.15 case shown in 
Fig. 2] We should note here that the AF moment extends 
outside the core, which is consistent with experiments. 

The LDOS obtained at the center of the vortex core 
is shown in Fig.^c). We can see that the LDOS shows 
neither the ZEP nor an explicit double-peak structure. 
To clarify the effect of the local AF moment on the LDOS, 




FIG. 3: The hole density n£ (a) and the LDOS (b) around 
the vortex core obtained by solving the BdG equations without 
AF order parameter. The doping rate is S = 0.11. 

we solve the BdG equations without AF order parameter 
for the same doping rate 6 — 0.11, as shown in Fig. [3] 
The LDOS obtained here (Fig. Efb)) is quite similar to 



that shown in Fig. Hfc), but at this time, a small peak 
can be found at E ~ 0. This peak is regarded as a 
trace of the ZEP. The only effect of the presence of AF 
order parameter is to annihilate the small peak at E ~ 0. 
We have confirmed that the spatial dependence of the 
amplitude of d- and s-wave components, and the hole 
density (Fig. Efa)) are quite similar to those shown in 
Fig. Ufa) and (b). 

From this comparison, we speculate as follows. The 
overall reduction of the spectral weight inside the vor- 
tex core already occurs without AF order parameters. 
Therefore this reduction is not due to the presence of the 
AF moment but due to the smallness of the core. For this 
doping rate (S = 0.11), the hole density inside the core 
decreases. This causes the appearance of the AF vortex 
core found in Fig. |2[b), but the effect of the AF core is 
subsidiary for the LDOS. 
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FIG. 4: The LDOS as a function of energy E/t and dis- 
tance r/a from the vortex core center, where a is the lattice 
constant. The range of the energy is —0.3 < E/t < 0.3. 

In order to clarify the origin of the reduced spectral 
weight of ZEP due to the smallness of the core, we study 
the spatial dependence of the LDOS for several hole dop- 
ings. Figure H shows the LDOS at S = 0.15,0.12 and 
0.11 as a function of energy E/t and distance r/a from 
the vortex core center, with a being the lattice constant. 
Here the distance dependence is plotted from the vortex 
core center (r = 0) along the (1,0) direction of the square 
lattice. First we look at the result for S = 0.15 shown 
in Fig. Ufa)- As r/a increases from 0, the ZEP first be- 
comes broader and smaller, then split into two peaks for 
r/a > 3. This kind of behavior of the LDOS, that is, the 
splitting of the ZEP away from the vortex core center, 
has been reported for s-wave vortex cores by Gygi and 
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Schluter[27j. They have shown that, in the s-wave case, 
the two peaks found away from the vortex core center 
correspond to the electronic bound states (core states) 
having higher angular momentum fi. The angular mo- 
mentum is quantized such that fi = m+ i where m is an 
integer. For small /i, the energy eigenvalues of the core 
states are known to have the approximate dispersion re- 
lation Efj_ ~ ujofi where luq is an energy quantum J^. In a 
semiclassical sense, the vortex core states have maximum 
amplitude at a distance r ~ fJ>/pF from the vortex core 
center. Of course in cZ-wave superconductors, the angular 
momentum is not conserved in a strict sense because ax- 
ial symmetry is broken. Quite recently, however, it has 
been pointed out that there is an alternative quantum 
number in d-wave superconductors which corresponds to 
the angular momentum)^]. Therefore the distance de- 
pendence of the splitting amplitude of ZEP shown in Fig. 
^indicates the dispersion relation of the core states: i.e., 
Efj, ~ tuofJ.- At S — 0.12 shown in Fig. |H[b), the split- 
ting of the ZEP exists even at r/a = 1, that is, the ZEP 
can be seen only at the vortex core center. Note that 
the AF moment is not induced for this case. Comparing 
the results in Fig. Ufa) and (b), we find that the slope 
of the dispersion relation becomes greater as the doping 
rate S decreases. This greater slope corresponds to larger 
luq or larger interlevel spacing of the core states, which 
is probably caused by the decrease of the core size. As 
the doping rate decreases further (6 = 0.11), the LDOS 
around E ~ is significantly reduced as shown in Fig. 

Finally let us comment on the effects of the next- 
nearest- neighbor hopping t'. The critical value of S where 
the AF core appears depends on the value of t' which re- 
inforces our viewpoint. For some values of t' , the AF core 
is realized even when the size of the vortex core is not so 
small. In this case, the spectral weight of the ZEP is not 
reduced owing to the large size of the vortex core, and 
there appears a large double peak structure in the LDOS 
caused by the presence of the AF moment [2ll |2^ | . This, 
however, does not explain the reduction of the spectral 
weight inside the core observed in YBCO or BSCCO. 

In summary, we have investigated the electronic states 
around vortex cores in high-T c superconductors based on 
the t-J model. The induced local AF moment and the 
LDOS without ZEP are found inside the vortex cores. 
However, the origin of the absence of ZEP is mainly due 
to the smallness of the vortex core sizes, which can be un- 
derstood within the conventional BdG theory with short 
coherence length. We find that the doping dependence is 
very important and thus it is quite interesting to study 
the overdoped region or low-T c superconductors experi- 
mentally, in which the vortex core size is large and the 
ZEP with the spatial dependence as shown in Fig. 4(a) 



will be observed. 
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